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ABSTRACT 

We present measurements of the UV galaxy luminosity function and the evolution of luminosity density 
from GALEX observations matched to the Sloan Digital Sky Survey (SDSS). We analyze galaxies in the 
Medium Imaging Survey overlapping the SDSS DR1 with a total coverage of 44 deg 2 . Using the combined 
GALEX+SDSS photometry, we compute photometric redshifts and study the LF in three redshift shells be- 
tween z = 0.07 and 0.25. The Schechter function fits indicate that the faint-end slope a is consistent with —1.1 
at all redshifts but the characteristic UV luminosity M* brightens by 0.2 mag from z = 0.07 to 0.25. In the 
lowest redshift bin, early and late type galaxies are studied separately and we confirm that red galaxies tend 
to be brighter and have a shallower slope a than blue ones. The derived luminosity densities are consistent 
with other GALEX results based on a local spectroscopic sample from 2dF and the evolution follows the trend 
reported by deeper studies. 

Subject headings: ultraviolet: galaxies — surveys — galaxies: luminosity function, evolution 



1. INTRODUCTION 

In the era of precision cosmology, the star formation history 
of the universe can be studied accurately as one can detect 
evolutionary effects in the observables on top of the global 
expansion. In particular, the restframe ultraviolet luminos- 
ity of galaxies has proven to yie ld a good handle on the star 
formation rate jKennicuttl 1 1 99 81) . A number of galaxy sur- 
veys have probed the history of star formation at different red- 
shifts. While most studies agree on a relatively rapid rise in 
the star formation rate (SFR) out to redshift of z ~ L s ignif- 
icant uncertainties remain even at lower redshifts (Lilly et aJJ 
| 1996| IConnollvetaIlll997l ICowie etaJ.lll999t IWilson et al. 
2002). The restframe UV continuum of local galaxies is not 
accessible from the grou nd. The balloon-b orne telescope of 
the FOCA experiment (Milliar ds all 19921) had been the best 
window onto the UV sky until last year, when the Galaxy Evo- 
lution Explorer (GALEX) satellite was successfully launched 
to orbit. This paper is one in the first series of luminosity 
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function papers on GALEX sources and focuses on galaxies 
at redshifts between z = 0.07 and 0.25. We use photometric 
redshifts to boost our sample size by a factor of 20 compared 
to spectroscopic data available. Throughout the paper, we as- 
sume a flat ACDM cosmology with = 0.3 and Hq = 70 km 
s" 1 Mpc" 1 . 

2. THE SAMPLE 

The GALEX telescope has two photometric bands at 1530A 
(FUV) and 2310A (NUV) and a 1.2 degree field of view. 
For a detaile d descriptio n of survey and per formance, see 
iMartin et all ( 120041) and lMorrissev et all ( 120041) in the present 
volume. Our sample consists of Medium Imaging Survey 
(MIS) fields overlappin g with the Data R elease One cov- 
erage of the SPSS i lDRl : Abazaiian et aj] 120031) . see also 
Seibert etaJ. (2004). The depth of the MIS fields is well 
matched to SDSS and this unique 7-band multicolor dataset 
provides a good basis for various statistical studies. We select 
the 57 fields with more than 1400 second exposure times and 
no objects with higher extinction than E(B-V) = 0.08. We 
use 36 arcmin radius circles in the center of the fields to en- 
sure uniform image quality. The intersection of the unique 
area of these MIS fields with the SDSS DR1 footprint is 43.9 
square degrees. Our catalog contains only objects that are 
classified as galaxies by the SDSS photometric pipeline based 
on their morphology. We use total magnitudes corrected for 
foreground extinction: SExtractor's MAG_AUTO for GALEX 
and model magnitudes from SDSS. For the limiting magni- 
tudes, we elect to choos e a safe m\\ m = 2 1 .5 cut in both bands 
to ensure completeness dXu et al . 200#. 

2.1. Photometric Redshifts 

Photometric redshifts are utilized to fully exploit the data 
set. We choose empirical photometric redshifts over template 
based estimates because currently the GALEX photometric 
system is only known to about 10% accurac y and SEP fit- 
ting is sensitive to zeropoint errors. Following Connoll y et all 
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FIG. 1. — Comparison of the spectroscopic and photometric redshifts for 
the MIS objects in our training set. The equation of the black dashed line is 
Zphot = 0. 

dl995al) . a third order polynomial formula was applied to map 
the GALEX NUV and SDSS u'g'r'i'z' magnitudes to photo- 
metric redshifts. Note that the FUV magnitude was excluded 
from the fitting formula as the NUV observations go deeper 
and not all galaxies have FUV measurements. This way there 
is only one redshift estimator that can be used for both FUV 
and NUV limited samples. We find that the empirical fit yields 
reliable redshift estimates out to redshift of 0.25. For the train- 
ing set of 6295 galaxies, the rms scatter is Azn„ s = 0.026 and 
there is about 2% outliers. This accuracy is about 15% better 
than SDSS alone using the same technique. A more detailed 
analysis and a photometri c redsh ift catalog will be published 
elsewhere (Budavari et aj] l2004cl) . We expect the uncertainty 
in the photometric redshifts to be a significant source of error 
in our statistical analysis, so we adopt a conservative nominal 
redshift error of a z = 0.03. Figure ^ compares the spectro- 
scopic and photometric redshifts as a function of redshift. 

One of the disadvantages of the empirical photometric red- 
shifts is that the method does not provide a direct measure- 
ment of the spectral types or K-corrections. To overcome this, 
we fit synthetic model spectra with full wavelength coverage 
from the UV to the IR to the SDSS photometry and pick the 
best fitting template for each galaxy. Our temp late set has 
10 inte rpolated spectra from Ell to Irr of Bruzua l"& CharloJ 
(2003). For each of these templates, the K-correction is cal- 
culated as a function of redshift. In addition to the galaxy 
templates, we also include a series of QSO spectra in an at- 
tempt to identify AGNs in the sample. Those objects which 
are best fitted with quasar templates, roughly 10%, are re- 
moved from the sample. Our photometric redshift catalog 
contains 190,489 MIS galaxies out of which 9,356 pass the 
area, magnitude, redshift and SED cuts in the NUV and 6, 174 
in FUV 

3. LUMINOSITY FUNCTION RESULTS 

There are several methods for ca lculating the luminosity 
function (iSchmidtl H%8t iLvnden-Belll 119711 iCholoniewskil 
119861 iSubbaRao etai.1 119961) . We use the VW method 
dSchmidtlll968r to calculate the LF in 0.1 magnitude wide 
bins, First, we derive the absolute magnitude using the dis- 
tance modulus and the K-correction, then the maximum red- 
shift where the object could be observed from. The LF is then 
calculated as <j>{M)dM = El/ViW), where V(z) = jd 3 (z) 
for a flat universe, f2 is the areal coverage and d(z) is the co- 
moving distance. 

To estimate the uncertainty in the LF, we create 50 Monte- 
Carlo (MC) realizations of the catalogs drawing the redshifts 
randomly from Gaussian distributions with means of the orig- 
inally estimated redshifts and widths of a z = 0.03. We coadd 



the MC realizations for a more robust estimate of the true LF. 
This method allows us to propagate the errors. The error- 
bars plotted in the figures are combinations of the variations 
among the MC realizations and the Poisson errors added in 
quadrature. 

3.1. Evolution with Redshift and Spectral Type 

To study the evolution of the UV LF as a function of red- 
shift, we split the sample into three redshift shells. These low, 
medium and high redshift subsamples have galaxies in the 
0.07-0.13, 0.13-0.19 and 0.19-0.25 intervals. FigureEJshows 
the FUV and NU V LFs for the thre e redshift slices along with 
their best fitting LSchechterl JT976h functions. The absolute 
magnitude range over the luminosity functions can be fitted is 
limited at the faint end by the lower redshift cutoff in the more 
distance shells and also at the bright end at M<-20, where the 
measurements depart from the Schechter function. The latter 
is due to residual contamination from QSO light that the SED 
fitting could not eliminate completely. The insets show the 
la, 2a and 3a confidence regions on the M* - a plane. As 
seen in Figure 12 there is a modest evolution in M* with red- 
shift in both bands but the leverage is not enough to constrain 
a to high accuracy at higher redshifts. In all cases, the slope 
is consistent with a = -1 . 1 . Table^lists the Schechter param- 
eters. 

We further divide the lowest redshift NUV and FUV lim- 
ited samples into two spectral classes based on the assigned 
SEDs. The early type galaxy class consists of objects with 
the five reddest templates and late type galaxies with the five 
bluer SEDs, which corresponds to a restframe color cut of 
(u'-r'\) = 1.7. This technique is expected to be more robust 
than the actual (V-r')o discriminator as all multicolor infor- 
mation is used. Figure|5]shows the FUV and NUV luminosity 
function for the early and late type galaxies. We find that M* 
is brighter for the red population by approximately 0.2 mag- 
nitudes in FUV and by 0.4 in NUV and a is shallower than 
for the blue population by 0.3, see Tabled The marginalized 
errors on both M* and a are large but the difference between 
the joint distribution of parameters is significant (see insets 
of Figure |3) When using the restframe colors to split galax- 
ies into red and blue, the plots exhibit the same features; the 
LFs at the faint end (M£— 18) are essentially indistinguishable 
and the bright end is also consistent with the template based 
results although the difference is slightly less pronounced as 
expected due to the larger scatter. The inverse concentration 
indices, simply the ratios of the radii containing 50% and 90% 
of the Petrosian r' fluxes, scatter significantly but on aver- 
age they are larger for galaxies in the blue class than the red 
ones, Cy ue = 0.46 and C~' d = 0.42. The red value is noticeably 
higher than that of t he elliptical SDSS galaxies, C" 1 w 0.35 
dStrateva et aJ.l20 01 ). which indicates that our redder class too 
contains spiral galaxies. 

3.2. Luminosity Density 

We derive the mean luminosity density (LD) by integrating 
the luminosity with the Schechter function, Pl = j 4>{L)LdL = 
<fi*L* T(a+2). In fact, the integral is calculated for not just the 
optimal fit but over the whole range of parameters. Weighting 
the results by the probability, w = exp(-\ 2 /2), is essentially 
the same as using Monte-Carlo realizations for estimating the 
errors, 5p\ = (p]} K - {pl) 2 k - Since the redshift range is limited, 
the errorbars on the higher redshift bins are rather large. Ta- 
ble^shows the LD measurements along with the statistical er- 
rors and estimates for the systematics due to cosmic variance 
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FIG. 2. — The GALEX FUV and NUV luminosity functions in three redshift bins {black open circle: 0.07 < z < 0.13, green solid circle: 0.13 < z < 0.19, 
magenta star: 0.19 < z < 0.25.) The top panel illustrates the 1/V ma x measurements along with the best fit Schechter functions. The confidence regions on M* 
and a are shown in the insets. The bottom panel shows the number of objects involved in the analysis for a particular selection. 




FIG. 3. — The luminosity functions differ in both passbands for the early and late type galaxies in the lowest redshift shell. The red galaxy population is brighter 
and has a shallower faint-end slope than the blue. 
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TABLE 1 

SCHECHTER PARAMETERS AND LUMINOSITY DENSITY 



Passband 


Redshift 


Type 


M* 


a 


\g<t>* [Mpc- 3 ] 




FUV 


0.07-0.13 
0.13-0.19 
0.19-0.25 


all 


-17.97±0.14 
-18.07±0.17 
-18.15 ±0.17 


-1.10±0.12 
-1.09 ±0.23 
-1.03 ±0.34 


-2.35 ±0.07 
-2.33 ±0.07 
-2.35 ±0.06 


25.51 ±0.02 ±0.06 
25.58 ±0.07 ±0.04 
25.61 ±0.12 ±0.03 




0.07-0.13 
0.07-0.13 


early 
late 


-17.98 ±0.20 
-17.74±0.19 


-0.80 ±0.19 
-1.12±0.17 


-2.74 ±0.08 
-2.47 ±0.09 


25.06 ±0.02 ±0.06 
25.31 ±0.04±0.06 


NUV 


0.07-0.13 
0.13-0.19 
0.19-0.25 


all 


-18.54±0.15 
-18.57±0.17 
-18.74±0.23 


-1.12±0.10 
-0.97 ± 0.20 
-0.99 ±0.35 


-2.38 ±0.07 
-2.33 ±0.07 
-2.39 ±0.08 


25.71 ±0.02 ±0.06 
25.74 ±0.05 ±0.04 
25.79 ±0.10 ±0.03 



0.07-0.13 late -1S.1 1 ±0.17 -1.09±0.14 -2.48±0.08 25.43 ± 0.03 ± 0.06 



that were derived similarly as in lWvder et all J2004). As seen 
in Figure [4] both pfuv and p^uv increase with redshift and 
are consistent with (1 +z) 3 as well as (1 + The errorbars 
in the figure are the combinations of the two sources of errors 
added in quadrature. They do not include errors from calibra- 
tion uncertainties of ^10% that may account for <5 lg pl = 0.04 
in both bands. 
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4. DISCUSSION 



Using 2dF red shif ts, the 
IWvderet all J2001 andlTre 



loca l GA LEX studies by 
2J2004) derived very con- 



sitent results with the present findings. In the corresponding 
redshift ranges, the Schechter fits are mostly within the 68% 
confidence regions with perhaps the one exception of M* in 
NUV, which seems to be brighter based on this photometric 
sample. The reason is that the different magnitude cuts in the 
UV yield slightly different galaxy populations and our sample 
has more redder galaxies, which makes M* brighter. 

The LP m easurements are also in good agreement with 
IWvder et 

ai]l|2004), 

and the observed trend is consistent with 
the GA LE X Deep Imaging Survey r esults bv lArnouts et all 
(2004) and Schi minovich et allJ2004l) probing the higher red- 
shift universe out to z = 1.2. The evolution, in turn, is also 
consistent with results by Wil son et"aH J2002I) . who find it to 
be proportional to (1 ' 7±L0 for galaxies out to z = 1.5 at 
restframe 2500A. 

Going from the observed UV luminosity function to a SFR 
function is complicated by dust and the fact that the FUV and 
NUV light trace stars forming on different timescales. We will 
address this problem in subsequent GALEX papers. 
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FIG. 4. — The luminosity density as a function of redshift in FUV (blue 
open circles) and NUV ( green filled circles) along with local GALEX results 
(squares) by [Wydgr et alj (20041) . The dotted and dashed lines correspond to 
(1 + z) 3 and (1 +z) L> , respectively, scaled to fit the lowest redshift bin of this 
study. 
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